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ABSTRACT 

The Atmospheric Imaging Assembly (AIA) telescope on board the Solar Dy¬ 
namics Observatory (SDO) provides coronal EUV imaging over a broader tem¬ 
perature sensitivity range than the previous generations of instruments (EUVI, 
EIT, and TRACE). Differential emission measure tomography (DEMT) of the 
solar corona based on AIA data is presented here for the first time. The main 
product of DEMT is the three-dimensional (3D) distribution of the local differen¬ 
tial emission measure (LDEM). While in previous studies, based on EIT or EUVI 
data, there were 3 available EUV bands, with a sensitivity range ~ 0.60 — 2.70 
MK, the present study is based on the 4 cooler AIA bands (aimed at studying the 
quiet sun), sensitive to the range ~ 0.55 — 3.75 MK. The AIA filters allow explo¬ 
ration of new parametric LDEM models. Since DEMT is better suited for lower 
activity periods, we use data from Carrington Rotation 2099, when the Sun was 
in its most quiescent state during the AIA mission. Also, we validate the para¬ 
metric LDEM inversion technique by applying it to standard bi-dimensional (2D) 
differential emission measure (DEM) analysis on sets of simultaneous AIA im¬ 
ages, and comparing the results with DEM curves obtained using other methods. 
Our study reveals a ubiquitous bimodal LDEM distribution in the quiet diffuse 
corona, which is stronger for denser regions. We argue that the nanoflare heating 
scenario is less likely to explain these results, and that alternative mechanisms, 
such as wave dissipation appear better supported by our results. 

Subject headings: Sun: corona - Sun: UV radiation - techniques: miscellaneous 
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Introduction 


Coronal heating is still a open research topic (Klimchuk 2015; [Cranmer 20091) - Different 
kinds of models have been proposed and these can be mainly classified in two broad classes: 
wave-driven models, in which wave dissipation transfers energy to heat the coronal plasma, 
and nanoflare models, in which the magnetic energy is released on a small spatial scale 
in an impulsive way by different possible physical mechanisms (magnetic reconnection, 
turbulence, instabilities, etc). To discriminate among the models, knowledge of the 
coronal differential emission measure (DEM) is important, as different models predict 
it in a different way. Standard 2D determinations of DEM are affected by line-of sight 
(LOS) integration ambiguities. The tomographic 3D reconstruction of the DEM, known 
as local-DEM or LDEM, removes this problem, and allows studying how the plasma is 
distributed in temperature within each tomographic grid voxel. 

Differential emission measure tomography (DEMT) allows reconstruction of the plasma 
parameters (electron density and electron temperature) in 3D, using temporal series of 
extreme ultraviolet (EUV) images. In a first step, a temporal series of narrow-band EUV 
images in a particular bandpass is used to determine the 3D distribution of the plasma 
emission in that channel, which we will refer to as filter band emissivity (FBE). The 3D FBE 
is found for every band of the EUV telescope independently. In a second step, the FBEs in 
each tomographic grid voxel are used to determine the LDEM, which describes the thermal 
distribution of the plasma within the voxel. Finally, moments of the LDEM are taken in 
each voxel to determine the mean electron density and the mean electron temperature. 
Previously published DEMT works are based on images from the Extreme Ultraviolet 
Imager (EUVI), on board the Solar TErrestrial RElations Observatory (STEREO) mission, 
a telescope that has 3 coro nal ban ds covering a temperature sensitivity range between 0.6 


and 2.7 MK. 


Vasquez et al. 


(120091 1 produced the first observational 3D analysis of coronal 
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prominence cavities, while IVasquez et al.l (120101.1201 lh analyzed the global st ructure of the 


solar corona during the minim um of activity between solar cycles 23 and 24. 


(12012(1 and 


Nuevo et ah 


Huang et ah 


(120131 ) combined DEMT with potential field source surface (PFSS) 


models of the global magnetic held to study the temperature structure of the quiet diffuse 
corona. Their studies revealed the ubiquitous presence of magnetic loops with downward 
gradients of temperature, and showed that they are present only during the solar minimum 


epoch. DEMT h as been also succesfully used as a y alidation tool for MH 


the solar corona ( Jin et ah 


2012 


Evans et ah 


2012) and the solar wind (lOran et ah 


3 modeling of th e 


2015|) 


The Atmospheric Imaging Assembly (AIA) telescope on board the Solar Dynamics 
Observatory (SDO) provides 6 coronal bands covering a temperature sensitivity range 
between 0.5 and 15 MK. The increased number of filters and their broader temperature 
range covered compared to EUVI, allows exploration of new models for the LDEM (or 
standard 2D DEM), providing more accurate estimates of the plasma parameters. As 
DEMT focuses on the slowly evolving quiet Sun, the AIA filters more appropriate for DEMT 
studies are those of 171, 193, 211, and 335 A, with maximum sensitivity temperatures in 
the range 0.9 to 2.5 MK. The 94 and 131 A bands have been primarily designed to record 
very hot flaring plasma, and thus exhibit a much smaller signal-to-noise (S/N) ratio in the 
quiet diffuse corona, which leads to more noisy FBE reconstructions. In the case of the 


94 A filter, its temverature response function 

TRF) includes a peak at about 7 MK and 

another one at about 1.25 MK ( 

Boerner et ah 

2012; 

Lemen et ah 

2012 

) but its sensitivity 


at the lower temperatures is poorly known (lAschwanden et ah 20131). These bands could 
be used as well, but this is deferred to a future publication. Finally, the 304 A emission, 
that corresponds to He II and Si XI, is optically thick on the disk and hence not useful for 
quantitative tomography (see Frazin et ah 2009). 


DEM inversion problems are difficult because of their ill-posed nature. When the 
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DEM is inverted using spectra as constraint, such as those provided for example by the 
Extreme-ultraviolet Imaging Spectrograph (EIS) on board Hinode, Monte Carlo Markov 
Chain (MCMC) approaches (Kashyap & Drake, 1998), or regularized inversion techniques 
(Hannah & Kontar, 2012) may be applied, with no need to prescribe a parametric functional 
form for the DEM. These methods have too many free parameters to apply to data produced 
by narrow band filters. The capa bilities of MCMC inversion techniques were recently 


investigated by 


Testa et al. 


( 2012 ) using synthetic EIS and AIA images computed from a 


3D radiative MHD simulation. They found that the inversion method is able to reproduce 
some global properties of the thermal distribution when using EIS images. The inversion is 
less accurate when using AIA images due to the more limited number of constra ints and the 


broad nature of the TRF of its filters. Similar conclusions were found by 


Del Zanna 


(12Q13|), 


who explored the limitations of applying MCMC methods to real AIA images. Also, using 
synthetic spectra computed from assumed DEM distributions, Landi et al. (2012) studied 
the ability of MCMC methods to reconstruct nearly isothermal plasmas, and found them 
unable to resolve them from multithermal plasmas with width smaller than AlogT = 0.05. 


In the case of narrow band images, parametrization of the DEM is a useful approach, as 
it provides simple, computationally efficient solutions. Parametric DEM studies, related to 
active regions, have been recently carried out by Aschwanden & Boerner (2012), Plowman 
et al. (2013), and Del Zanna et al. (2013). Parametric inversion has been implemented in 
DEMT to find the LDEM using the FBE values as constraints (Vasquez et al. 2009; Frazin 
et al. 2009). In this case, the values of the free parameters of the model are found as to 
minimize a functional that measures the mean square difference between the tomographic 
and synthesized FBEs. This is performed independently for each voxel of the tomographic 
grid. In previous DEMT works the LDEM was modeled by a single normal (Gaussian) 
function constrained by the FBEs provided by the 3 bands of the EUVI telescopes. In 
this work, a variety of parametric models for the LDEM is explored to best reproduce the 
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tomographic FBEs of the 4 AIA bands selected for DEMT. 

The article is organized as follows: In section 2, we describe the DEMT technique and 
the LDEM parametrizations that are tested. Sections 3.1 and 3.2 include a study of the 
LDEM inversion from simulated data and a comparison of DEMT results based on data 
from EUVI and AIA, respectively. Section 3.3 analyses the DEMT results using different 
parametric models. Section 3.4 shows an analysis of the bimodal properties of the LDEM 
in different coronal structures. Section 3.5 is 2D DEM validation study of the parametric 
inversion and the bimodal model. In Section 4 a discussion and summary of the results of 
this work is presented. 


2. Methods 

2.1. The DEMT technique 

DEMT consists of two steps. In a first step, a time series of EUV images is used to 
tomographically determine the 3D distribution the plasma emission in a given band (the 
FBE), formally defined in Equation (1) below. To that end the corona is discretized in 
a spherical computational grid. This first task is independently performed for each EUV 
band. In a second step, the FBE values of all bands are used as constraints to the inversion 
of the local temperature distribution (the LDEM) of the coronal plasma. This second task is 
independently performed for each voxel of the tomographic grid. From the resulting LDEM 
distributions, 3D maps of the coronal electron density and temperature can be derived. 
This section provides a brief description of the key concepts of DEMT that are needed in 
this work. For a detailed explanation the reader should consult Frazin et al. (2009). 

The AIA instrument images the off-limb corona up to about 1.25 i? Q . To perform 
EUV tomography, the inner corona volume in the height range 1.00-1.25 Rq is discretized 
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on a 25x90x180 (radial x latitudinal x longitudinal) spherical grid. Due to optical 
depth issues and EUV S/N levels, the DEMT results are reliable in the height range from 
1.03 to 1.20 Rq. For each band k, a series of images covering a full solar rotation is 
tomographically in verted to f ind, for each tomographic cell i, the value of the FBE Q , 


defined by fjFrazin et a! 


20091) 


cf^yWtAMA), (i) 

where r/j(A) is the plasma emissivity volumetrically averaged over the i-tli tomographic cell 
and (pk( A) is the normalized bandpass function of the band k. Due to temporal changes in 
the corona, tomographic reconstructions exhibit artifacts such as smearings and negative 
values of the reconstructed FBEs, or zero when the solution is constrained to positive 
values. These are called zero-density artifacts (ZDAs, see Frazin et al. 2009). Synthetic 
images based on the tomographic reconstruction can be produced by numerically computing 
line of sight (LOS) integrals of the 3D FBE. As an example, the left panel in Figure [D 
shows an image taken with the 193 A AIA filter (left) and the corresponding synthetic 
image computed from the tomographic FBE (right). The right panel shows the frequency 
histogram of the ratio of synthetic to observed intensity ratio for every corresponding pair 
of pixels in the two images. The relative difference between the synthetic and observed 
values in each pixel is below 10, 20 and 30%, for 37, 63 and 79% of the pixels, respectively, 
and virtually the same statistics hold when off-limb or on-disk pixels are considered 
separately. A similar level of agreement is found for all AIA bands. The tomographic model 
provides then a quite detailed reliable description of the average global corona during the 
reconstructed period. 

Within any given tomographic voxel, the LDEM £j(T) is a measure of the thermal 
distribution within the i th voxel. Following Frazin et al. (2009), the FBE and the LDEM 
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Fig. 1.— Left panel: image taken by the 193 A band of the AIA instrument (left), and cor¬ 


responding synthetic image computed by LOS-integration of the tomographic FBE (right). 
Right panel: frequency histogram of the synthetic to observed intensity ratio for every cor¬ 
responding pair of pixels in the two images. 

are related by 



( 2 ) 


where the temperature response function (TRF) Qk(T) is defined as 



(3) 


in which N e0 is a reference electron density and a 0 is a reference set of abundance used 
to computed the plasma emissivity from an optically thin plasma emission model such as 
CHIANTI (Dere et al. 1997). 

Due to the limited number of available bands, K, the inversion of LDEM function is 
under determined. The solution is implemented by modeling the LDEM with a family of 
functions £i(T) = F (T, af), depending on a vector of parameters c^, with L ~ K elements 
In each tomographic cell % the problem consists of finding the values of the parameters ccj 
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that allow to best reproduce the K tomographically reconstructed values of FBE in that 
cell. To do so an objective function is defined, that measures the quadratic differences 
between tomographically determined FBEs and the FBEs synthesized from the modeled 
LDEM, i.e. 


K 


<H ai ) = j2 c i 


-2 

k 


k =1 


Qk, ~ / dTQ k (T) F (T, cxi) 

^ ^"min 


(4) 


where K normalization coefficients C k , each given by the median value of Q k l in the whole 
reconstructed coronal volume, are included to make the objective function dimensionless. 
The minimization of the objective function is numerically implemented using the conjugate 
gradient method (?). 


At each cell i the following score is then computed, as a measure of the degree of 
success of the LDEM in reproducing the tomographically reconstructed FBEs 


Ri = (i/A')E|i-C„/c'S m |, 

k= 1 

where Ci,tom Y Cy'yn are the tomographic and synthetic FBEs. A perfect fit implies R 
and the higher the score the poorer the fit. 


(5) 

0, 


Once the LDEM is found at each tomographic cell i, we can derive 3D maps of the 
plasma parameters, by taking moments of the LDEM, specifically 


N 2 ■ = 

e,i 


T • = N~ 2 

- L m,i iV e ,2 


dT &(T), 


dT £i(T) T , 


( 6 ) 

(7) 


that are the square electron density and the mean electron temperature, respectively. 


The same formalism can also be used to perform a parametric 2D DEM analysis of 
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EUV image sets. The relationship between the intensity of band k in pixel i, and the 
DEM along the LOS associated with that pixel, ^(T) is, similarly to Equation ((2]), 

/ 7m ax 

dT Q k (T) ipi(T ), (8) 

min 

In the analysis below the DEM is described by the same parametric models used for the 
LDEM, and the optimal value of the parameter vector for each pixel is found in a similar 
fashion. In this case, the zeroth and first moments of the DEM distribution tpi(T) provide 
the emission measure (EM) and the mean electron temperature along the LOS associated 
with pixel i, respectively. In Section 13751 a parametric 2D DEM analysis of an active region 
(AR) and the diffuse region around it is developed, and the results are compared to those 
obtained by other authors. 


2.2. Temperature Response Functions 


Figure [2] shows the temperature response function (TRF) of the 4 coronal bands of 
the AIA instrument used in this work, as well as those of the 3 coronal bands of the 
EUVI. We used the latest effective areas released by the AIA team in the Solar Soft 
package, which includes a broad sensitivity plateau, roughly between 400 and 900 A, that 
affects the low temperature sensitivity of all bands. Calculations were performed with 


version 7.1 of the CHIANTI atomic data base ( Land i et a. 


set sun_coronal_f eldraan_1992_ext. abund (iFeldman et al 


2013b using the abundance 


1992; 


Landi et al. 


20020 and 


the ionization equilibrium calculations set chianti. ioneq. An important change in the 
spectrum introduced in CHIANTI v 7.1 is found in the wavelength range below 170 A, 
implying a more accurate calculation of the TRF of the 171 A band. With these settings 


we revisited the calcu 
previous calculations,! 


ation of the respon se function o: 


Lemen et al. 


2012 


Nuevo et al. 


’ all instruments. When compared to 


20121) the most important change is 


found in the very low temperature sensitivity of the AIA 335 A band, below the coronal 
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temperatures considered in this work. The most noticeable change at coronal temperatures 
is found in the TRF of the 171 band of both instruments, becoming narrower and thus less 
sensitive to temperatures below ~ 0.5 MK. All the TRFs are mainly characterized by a 
maximum sensitivity temperature (To)- Table Q] summarizes the values of To, as well as the 
full width half maximum (FWHMo). This information serves as a simple characterization 
of the sensitivity temperature range of each telescope band, and will be useful for setting 
up different aspects of the parametric LDEM inversion, as we describe in the next Section. 
Note that, in the case of the 335 A band, there is a secondary high sensitivity range that 
overlaps that of the 171 A. 


AIA & EUVI TRF [arbitrary unit] 



Te [MK] 


Fig. 2.— Temperature response functions of AIA (solid) and EUVI (dashed). 
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2.3. Filter Sets and LDEM Parametric Models 

In this section we describe and discuss the parametric models used to determine the 
LDEM. 


Single normal model (named “Nl r hereafter): 
«r) = yUV(T;[T 0 ,<7 T ]) = 


■ exp 


1 ( T-T 0 

2 l cr T 


(9) 


a/2tT0T 

In this model, the free parameters are the centroid To, the standard deviation ax and 
the amplitude A. This model will be applied when using the 3 AIA bands: 171, 193 
and 211 A a set of Liters we will refer to as “AIA-3” hereafter. This is the same 
parametric model used for all previous DEMT work based on data provided by the 
three bands of the EUVI instrument (e.g., Vasquez et al., 2011). The temperature 
range considered for the inversion, indicated as the limits in the integral of Equation 
([2]), is determined with the following criteria. The minimum temperature is taken to 
be roughly T min ps T 0 — FWHM 0 /2, where T 0 and FWHM 0 are the values of Table 
[U for the coolest band that is considered, or 171 A for AIA-3. Similarly, we take 
T max ~ T 0 + FWHM 0 /2, where T 0 and FWHM 0 are the values of Table Q]of the hottest 
band that is considered, or 211 A for AIA-3. In the case of both the AIA-3 and the 
EUVI sets this criteria leads to consider the temperature range [0.5, 2.3] MK. 

Double normal model (named “N2” hereafter): 


£(T) = A, A f(T- [To^an]) + A 2 A 7(T; [T 0>2 ,a T2 ]) , 


( 10 ) 


being a superposition of two normal functions. The first one has 3 free parameters (as 
for the N1 model). The second one has 2 free parameters (amplitude and centroid), 
with its standard deviation set at <jt 2 ~ 0.25 MK (which is half the temperature 
step between the peak response of the consecutive AIA bands 211 and 335 A). This 
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model has then a total of 5 free parameters, and will be applied when using the 4 
AIA bands: 171, 193, 211 and 335 A, a set of filters named ”AIA-4” hereafter. In this 
case, even if the number of free parameters is larger than the number of constraints, 
the approach is justified as it is posed as a global optimization problem. To select 
the temperature range to be considered for the LDEM inversion we apply the same 
criteria applied to Nl, which leads to consider the range [0.5, 3.8] MK. The strategy 
of this N2 model is to allow the first fully-free normal distribution to mainly describe 
the plasma below about 2 MK (detected by the 171, 193 and 211 A bands, as well as 
by the low temperature end of the TRF of the 335 A band), while the second normal 
aims at mainly describing the plasma with temperatures around 2.5 MK (detected by 
the main sensitivity peak of the 335 A band). Accordingly, the standard deviation 
value set for the second normal function is of the order of the temperature difference 
that the “neighbouring” 211 and 335 A bands can resolve. 

• Quadruple normal model (named “N4” hereafter): 

4 

t(T) = Y.AjM(T-,[T { , ( 11 ) 

3 =1 

being a superposition of four normal functions with their centroids and widths set 
at prescribed values, with 4 free parameters that are the amplitudes of each normal 
function. The fixed centroids are set to the T 0 values of the four TRFs (see table [1]) 
and the fixed widths are set equal to half the temperature difference between the To 
values of consecutive bands. The strategy of this parametrization, applied when using 
AIA-4, is for each normal function to predict the emissivity of the respective band. 

• Other unimodal models: For AIA-4 we experimented several alternate unimodal 
models, specifically: 


a) The Nl model. 
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b) The “top hat” (TH) model, 


£(T) = TH (T; [A, T 0 , a T , a]) = 


A 


■ exp 


1 f T-T 0 

2 V °t 


2 a 


( 12 ) 


\ZTkot 

characterized by 4 free parameters: an amplitude A, a centroid T 0 , a width 
a t, and a parameter a controling the gradient of the sensitivity drop at both 
the cooler and the hotter temperature ends of the model. Note that for a — 1 
the model becomes normal, while for larger values of a the shape of the model 
resembles that of a top hat function. 

c) The asymmetric “top-hat” (ATH) model, 

£(■ T) = ATH (T; [A, T 0 , a T , a, p]) = TH (T; [A, T 0 , a T , a}) T p , (13) 

which multiplies the TH model by a power law, adding a 5 th free (signed) 
parameter p that allows to achieve asymmetrical distributions. 


Figure [[2] shows examples for all the models described above. As mentioned in 
the int roduct ion, similar models have been used by recent parametric DEM studies. 


Aschwanden & Boerner 12 0111 ) performed DEM analysis along bright loops in ARs based 


on AIA images, modeling the DEM with a combination of different number of normal 
distributions in order to reproduce the emission observed in the six bands of AIA. Del 
Zanna (2013) also used AIA data to study the DEM in the core of an AR, modeling the 
DEM in each pixel with ~ 10 normal distributions spaced across the temperature range 
from 0.5 to 4.0 MK. 
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3. Results 

3.1. Inversion of Simulated Data 

With the aim of performing a controlled study of the LDEM inversion, and 
understanding the capabilities and limitations of the parametric models, we implemented 
an ensemble of assumed LDEM curves and derived synthetic FBE data. To generate the 
synthetic data we adopt for the LDEM a “top-hat” model with a very steep gradient at 
both ends of the distribtion. The free parameters of the model are set so that the resulting 
distribution is roughly constant between a fixed minimum temperature ~ 0.8 MK (a coronal 
temperature a bit lower than the maximum temperature sensitivity of the 171 A band), 
and a maximum temperature that we change for each model. The ensemble is generated by 
changing the maximum temperature of the simulation, that was set at a minimum value 
of ~ 1.2 MK, and incremented in steps of 0.01 MK up to maximum value of ~ 10 MK. 
The simulation with the lowest maximum temperature has all the plasma in the sensitivity 
range of the 171 A band, while the successive simulations add plasma at progressively larger 
temperatures, in the range of sensitivity of all the bands of the AIA-4 filter set, and beyond. 

For every assumed LDEM of the ensemble we calculated its moments using Equations 
©-(0). Using the TRFs showed in section EOl we calculated synthetic values of the FBE 
for all bands using Equation (|2]). We then used this synthetic data to perform the LDEM 
inversion, as described in section I2TT1 using the parametrizations N1 and N2 explained in 
section 12.31 Finally, we calculated the moments of the inverted LDEM. We are now in 
a position to compare, for all simulations in the ensemble, the moments of the assumed 
LDEM and those of the inverted LDEM. Note that the way the simulations are designed, 
the mean temperature T m of the simulations monotonically increases with the density N e , 
which is important to bear in mind when analyzing the results below. 
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Using the AIA-3 data set and the N1 model for the LDEM (AIA-3/N1, hereafter), the 
top panels of Figure [3] show scatter plots of the moments of the assumed LDEM versus 
those of the inverted LDEM. The bottom panels show the same comparisons when using 
the AIA-4 data set and the N2 model for the LDEM (AIA-4/N2, hereafter). 

Ne [x 10 s cm -3 ] Tm [MK] 



0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 
assumed 


Ne [x 10 8 cm 3 ] 



Tm [MK] 



Fig. 3.— Scatter plots of the assumed LDEM moments versus the inverted LDEM moments 
for AIA-3/N1 (first row) and AIA-4/N2 (second row). 

For AIA-3/N1, the simulations with T m <1.7 MK show a very high correlation between 
the moments of the simulated and the inverted LDEM, while for simulations with T m >1.8 
MK the moments of the inverted LDEM are approximately constant. The reason for this is 
that for top-hat LDEM simulations with T m >1.8 MK the maximum temperature becomes 
Xmax > 2.8MK, so that the plasma that is added by models with increasing T m is well 
beyond the upper limit of the high sensitivity range of the 211 A band, which, from Table 
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1 is roughly To + FWHMo/2 ~ 2.3 MK. 

Similarly, for AIA4/N2 (second row of panels) there is a good correlation between the 
moments of the simulated and the inverted LDEM for simulations with T m < 2.7 MK, while 
for simulations with T m > 2.8 MK the moments of the inverted LDEM are approximately 
constant. The reason for this is that for top-hat LDEM simulations with T m > 2.8 MK 
the maximum temperature becomes Tmax > 4.8MK, so that the plasma that is added by 
models with increasing T m is well beyond the upper limit of the high sensitivity range of 
the 335 A band, which, from Tabic 1 is roughly T 0 + FWHM 0 /2 ~ 3.75 MK. 

As previous DEMT works have been based on the EUVI instrument, we are interested 
in testing how the results of inverting the LDEM based on data from that instrument 
compares with those of inverting the LDEM based on data from AIA-3. Towards that 
end we have also performed a similar analysis of inversion of simulated data, as described 
above, but using the TRFs of the EUVI-B instrument and the N1 model for the LDEM 
(EUVI-B/N1, hereafter). For every simulation of the ensemble, we compared the moments 
of the LDEM inverted from EUVI-B data against those of the LDEM inverted from AIA-3 
data, using the N1 model in both cases. The inferred electron density differ in less than 3% 
in all cases, while the mean temperature differ in less than 5%. 


3.2. DEMT Results with Three Bands 

We performed the DEMT analysis of the Carrington Rotation (CR) 2099 (2010, 13 
July through 9 august), during the early rising phase of solar cycle 24. This period is one 
of the earliest ones for which AIA data is available for the full rotation, and the most quiet 
period ever observed by that instrument, which is preferred for global DEMT analysis. 

As all previously published work on EUV tomography has been based on EUVI data, we 
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also performed reconstructions based on EUVI data for the same period, for comparison 
purposes. During CR-2099 the angular separation between SDO and STEREO was about 
~ 72°. As a consequence, differences observed between the results using AIA and EUVI 
are partly caused by solar dynamics. As the DEMT analysis is not well suited to study 
ARs, these were excluded from it. In the quantitative analyses shown below, tomographic 
voxels belonging to the open corona, the closed quiet sun (QS) corona, and the ARs where 
carefully separated. To this end, ARs were identified in the catalogue provided by the 
National Oceanic and Atmospheric Administration (NOAA) Space Weather Prediction 
CenteiQ, and verified that they correspond to the highest density values in the tomographic 
reconstruction. Indeed, we were able to indentify all ARs by finding the regions where the 
tomographic density value was above a suitable threshold value. A threshold was established 
at each height of the tomographic grid independently, since the density decreases rapidly 
with height. For example, at 1.075 R Q a good way to characterize ARs in the tomographic 
reconstruction is to look for tomographic density values above ~ 1.4 x 10 8 cm -3 . 

Figure |4] shows the DEMT results at a sample height of 1.075 R & using both 
instrumental sets (with 3 bands) and the N1 LDEM parametrization: EUVI-B/N1 and 
AIA-3/N1. Similar maps are also obtained at all 25 height bins of the tomographic 
grid. Using a PFSS extrapolation we overplot the boimdaries between magnetically 
open and closed regions as a thick black curve. The PFSS model is computed using the 
finite-difference iterative solver FDIPS by Toth et al. (2011) on a 150 x 180 x 360 spherical 
grid, covering 1.0 to 2.5 Rq. While the morphology of the corona is clearly more complex 
compared to that of solar minimum, there is an overall good agreement between the 
open/closed boundary of the PFSS model and the density and temperature structures of 
the DEMT analysis. There is an overall high consistency between the results of the PFSS 

1 www.swpc.noaa.gov/ftpraenu/warehouse/2010.html 



19 


and the tomographic models. Figure® shows that the location of the open/closed boundary 
of the PFSS model is characterized by a very high transverse gradient in both the electron 
density and the mean temperature maps derived from the DEMT analysis. This has been 

yzed in detail in DEMT studies of solar minimum rotations 
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In the mean temperature maps, black voxels indicate ZDAs, the regions of zero 
emissivity described in Section 12.11 In the same maps, white voxels indicate cells for which 
the parametric LDEM possesses a score R > 0.2, in these cases, we consider that the 
parametric LDEM performs poorly in reproducing the tomographic FBEs. This happens in 
locations where the set of FBE values is unusual, and we dub these anomalous emissivity 
voxels (AEVs). In the electron density maps, both ZDAs and AEVs are indicated as black. 
In the score R maps, ZDAs are indicated as black, while AEVs are not highlighted. 


Before comparing the similarities and differences between the DEMT results obtained 
with the different instrumental sets, we will comment on the quality of the results as 
measured by the score R defined in Equation (5). Figure® shows histograms for the values 
of R shown in the bottom panels of Figure ® for the open and closed region separately. 
The vast majority of the tomographic voxels possess a score R < 10~ 2 , meaning that the 
tomographic and synthetic band emissivities differ less than 1%, in average. 

As seen in Figure® the results obtained with EUVI/N1 and AIA-3/N1 are very 
consistent in the closed magnetic region. To evaluate differences quantitatively, Figure ® 
shows scatter plots and histograms comparing the DEMT results obtained with EUVI/N1 
and AIA-3/N1 in the closed region. There are small systematic differences, with a mean 
temperature and density increase of order ~ 5% for the AIA results respect to those of 
EUVI. The differences are due to differences in the TRFs of both instrumental sets, mainly 
between the AIA 211 A band and EUVI 284 A band TRFs. There are also non-systematic 














differences, reflected by the standard deviation of the histograms, which are to be expected 
due to solar dynamics effects as the time series of data used for each instrument are not 
synchronous. 


3.3. DEMT Results with Four Bands 

In this section we analyze the results of DEMT based on the AIA-4 data set, which 
adds the 335 A band to the AIA-3 set. The additional band has a maximum sensitivity at 
temperature 2.47 MK. 

First, the LDEM is modeled with the N1 parametrization. The results are displayed in 
the left column of Figure [7] which shows, from top to bottom, the Carrington maps of N e , 
T m and R at 1.075 Rq, and the frequency histograms of the R values at the same height, 
separated in the open and closed region. A comparison of the results for R in this column 
against those in the right column of Figure SJ and also those in Figure 0 shows that the N1 
model can not explain the emisivity in the 4 bands as succesfully as it does for 3 bands. 

It is the shape of the N1 model and the fact that it depends on only 3 parameters which 
makes it unsuccessful in coping with the wider temperature sensitivity range of the AIA-4 
data set. 

Next, we test the TH parametric model, which by adding a 4 th parameter has the 
flexibility of widening the range of temperatures over which the LDEM is near its maximum 
value, as compared to the N1 model, while still being a simple symmetric unimodal 
distribution. The middle column in Figure [7] shows the results using this model. Even if 
the number of free parameters has been increased by 1 compared to the N1 model, the 
histograms of R do not show improvement over the previous model. 

Next, we test the ATH model, which allows for both negative and positive temperature 
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gradients, depending on the sign of the exponent of the multiplicative power law, allocating 
more plasma at the low or high temperature end, as needed. The results of using this 
unimodal assymetric parametrization are shown in the right column of Figure [7] The 
histograms of R show a slight improvement over the N1 and TH models, with the smallest 
median and mean values for the score R , while still performing significantly worse than the 
N 1 model for the AIA-3 set. 

Figure [7] shows that, in order to explain the emissivities of the AIA-4 set none of 
the unimodal models, symmetric or asymmetric, performs as well as the N1 model does 
for the AIA-3 set (Figure [I]). Looking for better agreement, we attempt the multimodal 
distributions N2 and N4 described in Section 12731 The maps of N e , T m and R are shown in 
Figure [H and the histograms of R are shown in Figure |9j 

A qualitative comparison between the results of the N2 and N4 models (left and right 
panels of Figure [ 8 ]) shows that both models have similar values of N e and T m . Indeed, a 
quantitative comparison (not shown to save space) reveals that the median LDEM moments 
differ in less than 4%. Even so, the N2 model performs better than the N4 model, as it 
achieves the best score R. Indeed, a comparison of the histogram of the left panel of Figure 
[9] with the histograms of Figure [5l reveals that the score R of the N2 model when using the 
AIA-4 data set is as good as the one achieved by the N1 model when using the AIA-3 data 
set. Figure flOl shows the histograms of the ratio ^ s y nth )/^h om ) f or the 4 b anc i s k of AIA-4, 
for the different parametric models used for the LDEM, at 1.075 R Q . It is readily verified 
the much higher degree of suceess of the N2 model for each band independently. 

The high degree of success of the N2 parametrization of the LDEM to reproduce 
the FBEs observed with all 4 bands of the AIA-4 instrumental set is due to the assumed 
bimodality. Our results show then that the quiet solar corona local (3D) DEM is multimodal 
at the spatio-temporal resolution of the DEMT and the temperature resolution of the TRFs 
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of the AIA filters. 

The differences seen in Figures [3] and O between the results obtained with AIA-3/N1 
and with AIA-4/N2 data sets are mainly due to the inclusion of AIA 335 A band, sensitive 
to temperatures higher than the other three bands. As can be qualitatively seen in the right 
panels of Figure® and the left panels of Figure® the AIA-4/N2 results have systematically 
larger values of N e and T m respect to the results obtained with AIA-3/N1. The systematic 
increase in the LDEM moments is more significant in the closed region, where the hotter 
plasma is located. To evaluate differences quantitatively, Figure dU shows scatter plots 
comparing the plasma parameters obtained in the closed region with N2/AIA-4 and those 
obtained with Nl/AIA-3, as well as histograms of the respective ratios. When the 335 A 
band is included the LDEM inversion leads to an increase in the electron density and mean 
temperature of about 14 and 18%, respectively, in the closed region. A similar analysis, but 
considering the voxels in the open region (not included here to save space) indicates that 
the increase in the density and temperature are about 3 and 10%, respectively. 

So far we have compared global maps of the moments and score R that are obtained 
with all the different parametrizations. It is interesting to also compare all the parametric 
distributions. For two sample tomographic voxels, one belonging to a hotter and denser 
coronal region than the other, Figure [®] shows the LDEM obtained using the different 
parametric models, and compares their performances. 

We first highlight that the Nl/AIA-3 (dashed-blue) and N2/AIA-4 (solid-blue) results 
achieve R < 10~ 2 in both coronal regions. Among the rest of the models, the one that best 
performs for AIA-4 in both voxels is the ATH model (solid-green). The reason is that by 
adjusting a positive or negative temperature gradient, this model is able to allocate more 
plasma in the large temperature end for the hotter cell, and in the low temperature end 
for the cooler one. In the case of the N1 model for AIA-4 data, the single normal adjusts 
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its centroid to an intermediate temperature between the two modes of the N2 model, but 
it is of course unable to allocate plasma preferently in one temperature end or the other. 
Interestingly enough the N4 model does not perform so well. The reason for this is that this 
model can only adjust the 4 free amplitudes to allocate more or less plasma at each fixed 
centroid temperature, which affects neighboring bands as their respectives TRFs overlap. 
The reason the N2 model performs so well is that the three bands 171-193-211 A are well 
explained by the cooler component, while the sensitivity at the maximum sensitivity of 
the 335 A band is explained by the hotter component. The low-temperature sensitivity 
in the 335 A band (that overlaps that of the 171 A band) is also captured by the cooler 
component. 

Finally, it is very interesting to note that the electron density and mean temperature 
predicted by the N2 model are almost the same as the values predicted by the ATH model, 
over the whole corona. This is shown in detail in Figure [13l which compares the electron 
density and mean electron temperature obtained with both models. Even if the choice of a 
reasonably low R parametrization is not unique, they achieve very similar 3D maps of the 
electron density and mean temperature. In this regard, the AIA-4 set is able to reveal extra 
plasma not previously detected by AIA-3, EUVI, or EIT, which accounts for an extra 14% 
of coronal mass in the closed region, and a mean temperature 18% higher. 


3.4. Analysis of the Multimodal Corona 

In the previous section we established that the closed quiet corona is bimodal at the 
spatiotemporal resolution of the tomographic grid and the temperature resolution of the 
AIA narrow band TRFs. In this section, we quantify how strong the bimodality of a given 
LDEM distribution is, by calculating the relative plasma fraction (N e j/N e ) 2 for each normal 
component (j = 1,2) of the N2 model. Figure fl4l shows Carrington maps of both ( N e ^/N e ) 2 
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and (TVg^/iVe) 2 , at 1.075 R Q . It is readily seen that the cooler component is dominant in 
most of the magnetically closed corona, except in localized regions, and in coronal holes. 
The hotter component is significant in the magnetically closed corona, and dominates the 
compact regions where the cooler component is more modest. To quantify the statistical 
relevance of the hotter component, Figure H5l shows the frequency histograms of (IVe^/iVe) 2 
at the same height in the open and closed regions, separating the voxels within the QS 
(middle panel) from those in the ARs (bottom panel). It is readily seen that the bimodality 
is much more important in the closed region than in the open region. Also, within the 
closed region, the bimodality is more important in ARs than in the QS. To quantify this, 
the fraction of voxels for which the hotter component is the dominant one, i.e. for which 
(N e ^/N e ) 2 A 0.5, is 48, 6 and 1 % in the ARs, closed quiet corona, and open corona, 
respectively. 

Figure fl6l shows Carrington maps and histograms of the centroid T 0 ,i °f the cooler 
component and the centroid T 0j2 of the hotter one. The values of T 0 ,1 are significantly and 
sistematically smaller in the open region, consistent with sub-MK temperatures in coronal 
holes reported by many studies (Vasquez et al. 2012; Vasquez et al. 2011; Feng et ah 2009; 
Wilhelm, et ah 1998). On other hand, the values of the centroid T 0j 2 distribute always 
around the maximum temperature sensitivity of the 335 A band (see table CO) throughout 
the closed corona. Within the open corona the centroid T 0i2 tends to be systematically 
smaller, yet statistically close to the maximum temperature sensitivity of the 335 A band. 
As commented above in relation to Figure d5j the hotter component of the model N2 is 
much weaker in the open region. This is due to the colder temperatures that dominate the 
open corona and the correspondingly much weaker 335 A emissivity. Still, the model N2 is 
the most successful one also in the open corona, as it allows to make the hotter component 
as weak as needed, yet still present. Unimodal models can mimic this by increasing their 
width, at the cost of diminishing the synthetic-to-observed agreement for the colder bands, 
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then increasing the value of the score R. 

In order to gain understanding of which locations exhibit more or less bimodality, 
we correlated the results of the AIA-4 and AIA-3 studies in the quiet corona. The left 
panel in Figure [171 shows a scatter plot of (N e ^/N e ) 2 obtained with AIA-4/N2 versus the 
electron density value N e obtained with AIA-3/N1. The orange dots indicate voxels of 
the magnetically closed region, while blue dots indicate voxels of the magnetically open 
region. The right panel shows a scatter plot of ( N e ^/N e ) 2 obtained with AIA-4/N2 versus 
the electron mean temperature value T m obtained with AIA-3/N1. These plots show that, 
in the closed quiet corona, the bimodality is stronger for denser and hotter regions, it 
also shows that when the cooler component is denser the same happens with the hotter 
component. As we will see in the discussion section this fact places an important constraint 
on heating models of the quiet corona. Finally, in the open region the same trends are 
observed and, in addition, the bimodality is on the average much less important than in the 
closed region, as already shown in Figure fl5l 

For the analysis shown in this section we chose 1.075 R Q as an intermediate sample 
height. Similar results were found at all tomographic heights. We investigated the 
dependence with height of both N e)2 {r ) and A r e l (r) within the quiet region of the streamer 
belt, similar to Vasquez et al. (2011). The scale heights of both quantities differ by 
only about 12%, so that the ratio (A^/A^q) 2 is roughly constant with height, showing a 
variability of less than 10% over the range 1.025 to 1.20 R & . 

Spectral observations of the corona allow studying the temperature of the emitting 
plasma by means of the emission measure (EM) technique. If the number of available 
spectral lines is large enough, the EM technique is a powerful tool to assess whether a 
plasma is isothermal or multithermal. Systematic studies of the temperature structure of 
the solar corona using this technique have been recently reviewed by Feldman & Landi 
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(2008). They conclude that the existing EM analysis of the corona support a scenario in 
which the coronal plasma is a combination of a few nearly-isothermal populations, with 
well defined temperatures. They find coronal holes to be characterized by temperatures 
around log(T) = 5.95, quiet sun regions by temperatures around log(T) = 6.15, and ARs 
appear characterized by two distinct populations with temperatures log(T) = 6.15 and 
log(T) = 6.5. While these works did not show bi-modal distributions, the temperatures 
reported therein are consistent with our results. According to the histograms in the 
bottom panels of Figure [HI the closed quiet sun is characterized by two well defined 
temperatures given by the mean values log < T 01 >= 6.15 and log < T 0i2 >= 6.41. The 
mean temperature of the cooler component of the ubiquitous bimodal LDEM found by the 
DEMT study matches the characteristic temperature of the quiet sun reported by the EM 
studies. The mean temperature of the hotter component of the LDEM matches that of the 
hotter component of the EM studies in ARs. 

Those EM studies where there is only a warm component (close to To,i) may not 
have included enough hot lines, or may have been centered on regions where the hotter 
component (close to T 02 ) was smaller. A weaker signal of the hotter component may also 
be the result of the LOS-integration encompassing more regions with the cooler component 
being dominant over the hotter component, giving a cumulative larger importance to the 
former than to the latter, which then does not provide enough counts to be observed. The 
case may also be that the reported EM studies of quiet sun regions correspond to epochs 
when the hotter component found in this DEMT study (that corresponds to the early 
rising phase of solar cycle 24) was particularly weak. In this regard, Landi & Testa (2014) 
find that overall the thermal structure of streamers does not change much over the solar 
cycle, although with a few spectroscopic measurements. LOS-integrated spectral studies 
may also be consistent with more than one characteristic temperature (Ralchenko et ah, 
2007). The question rises in that case of whether such multithermality is an artifact of the 


LOS integration of different large-scale structures at very different places. In this context 
it should be highlighted that in the results of this paper the multithermality is a local 
property. 


3.5. 2D DEM analysis 

The technique for the parametric inversion of the LDEM used in this work is the same 
as has been used for all the previously published DEMT work (for example Frazin et al. 
2009, Vasquez et al. 2010, 2011, and Nuevo et al. 2013). As explained in Section [2TTT the 
same technique can be used to find the standard 2D DEM from simultaneous EUV image 
sets. With the aim of cross-validating our parametric inversion technique, we use it to 
perform a standard DEM analysis of an AR and the diffuse region around it. We used the 
AIA-4 images taken at 12:00 UT on 2010, July, 25, and found the DEM at each pixel using 
the N2 model. By taking the zeroth moment of each DEM we find the emission measure 
(EM [cm” 5 ]) of each pixel. 

The left panel of Figure fl8l shows, in gray scale, an image of the EM of the analyzed 
region, which corresponds to AR NOAA 11089. The color boxes, of size 5x5 pixels, identify 
selected regions in the image. Using the same color code, the right panel shows the average 
DEM model obtained within each color box. This figure can be directly compared with a 
similar analysis shown in Figure 6 of Narukage et al. (2014), which performed a MCMC 
inversion. Their technique to determine the DEM, described in Winebarger et al. (2012), 
uses 7 spectral lines of the Extreme-ultraviolet Imaging Spectrograph (EIS) on board the 
Hinode (1-6 MK) and 2 broadband images of the X-Ray Telescope (XRT) on board Hinode. 
In their Figure 6, the authors show the DEM obtained in each selected region, which are 
indicated in their Figure 7. Note that we have chosen the color code of this work in a 
similar fashion to that of Figures 6 and 7 of Narukage et al. (2014), with the red/yellow 
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colors corresponding to the core of the AR core, dark/blue colors corresponding to the dark 
region around the AR, and green/light-blue colors corresponding to intermediate regions. 

There are several similarities between both works that we should highlight. In the 
comparison one should bear in mind that the specific ARs analyzed in both works is not 
the same. Firstly, the MCMC analysis of Narukage et al. (2014) shows in a bimodal DEM 
in virtually all regions, except for pixels in the deep core of the AR, where a third hotter 
component may arise. Secondly, in both analyses, as one moves from the dimmer regions 
surrounding the AR into the deep core of the AR, the area of both the cooler and the 
hotter component increases, the centroids of the two components gradually shift to larger 
temperatures, and the hotter component becomes dominant. In both works, the range of 
temperatures of the plasma described by the DEM exhibits similar values. For the MCMC 
analysis the centroid of the cooler component ranges from 1.2 to 1.5 MK, while the centroid 
of the hotter component ranges from 2.5 to 4.0 MK. In the parametric analysis of this 
work, the centroid of the cooler component ranges from 1.5 to 2.25 MK, while the centroid 
of the hotter component ranges from 2.5 to 4.0 MK. The absolute values of the DEM in 
units of [cm -5 K -1 ] (that we will omit in the rest of this paragraph) are also similar in both 
analyses, with comparable orders of magnitude. In the MCMC analysis, the amplitude of 
the cooler component is typically within the range 1 — 4 x 10 21 , while the amplitude of the 
hotter component ranges from 0.5 — 10 x 10 21 . The parametric analysis of this work, the 
amplitude of the cooler component ranges between 0.3 — 6 x 10 21 , while the amplitude of 
the hotter component ranges between 0.2 — 70 x 10 21 . In the core of the AR, the DEMs by 
Narukage et al. (2014) predict the intensity of the key EIS spectral lines (Winebarger et al., 
2012) within a 30% precision. That level of agreement is comparable to the requirement 
R < 0.2 we have set for the voxels that are included in our tomographic analysis. 
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4. Discussion and Conclusions 

In this work the DEMT technique is applied for the first time to 4 coronal bands 
instead of the usual 3 bands used in all previous works. The 4 bands used are those of 
the SDO/AIA instrument that are more suitable for studying the quiet corona, namely 
the bands of 171, 193, 211, and 335 A. The analysis of simulated data (Section l3.ip shows 
how the inverted LDEM is able to succesfully reproduce the electron density and mean 
temperature of an assumed LDEM (derived from its zeroth and first order moments), even 
when both LDEM parametric models are not the same. It also clearly illustrates how the 
EUV bands are incapable of detecting plasma at temperatures outside their sensitivity 
range, and thus how the LDEM aims at describing the plasma only in that sensitivity range. 

This work is the first to apply DEMT to AIA data and, therefore, a comparison of 
results against those based on EUVI data is included. The analysis indicates that there is 
a high degree of consistency between the DEMT results obtained with AIA-3 and those 
obtained with EUVI. DEMT quantitative results based on AIA-3 data can then be directly 
compared with previous published DEMT reconstructions based on EUVI data. 

Using the AIA-3 data set, the single normal (Nl) parametric model of the LDEM 
(3 free parameters) is able to succesfully reproduce the tomographic emissivities of all 
bands (171, 193, and 211 A), but it can not expain the AIA-4 data (adding the 335 A 
band). Other symmetric and asymmetric unimodal LDEM parameterizations were tested 
for AIA-4. In all cases the tomographic versus synthetic agreement is not nearly as good as 
when using 3 bands, despite having 4 or 5 free parameters. 

Multimodal LDEM were modeled using a variety of combinations of partially 
conditioned normal functions, an approach similar to Aschwanden & Boerner (2011). The 
model that consistently achieves the best goodness-of-fit score R is the bi-normal (N2) 
model, which is a superposition of two normal distributions with distinct cool and hot 
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components at (Top) ~ 1.4 MK and (T 0) 2 ) ~ 2.6 MK. Our study shows that the quiet 
corona LDEM is multimodal at the spatial resolution of the tomographic grid, which is 
0.01 R q x 2° x 2°, or about (7 x 10 3 km) x [(2.44 x 10 4 km) 2 ] for a representative voxel at a 
height of 0.1 Rq above the phostosphere at the equator. 

The analysis in Section 13.31 for AIA-4 shows that any choice of parametrization that 
reasonably reproduces the emissivity in all bands leads to very similar estimates of both 
density and mean temperature (see Figure fl2l) . which is often the reason to perform DEMT. 
Using the AIA-4 set it is found that the 335 A band reveals plasma previously not seen by 
instruments of 3 bands (EUVI), which provides for an increase of about 14% of the mean 
electron density in the closed quiet corona when compared to results with AIA-3, as well as 
for rise of about 18% of the estimated mean temperature. 

The bimodality is found to be stronger for denser and hotter regions. The cooler 
component of the bimodal LDEM is dominant in coronal holes and most of the magnetically 
closed corona, except in compact active regions where the cooler component is more 
modest. The fraction of voxels for which the hotter component is the dominant one, i.e. 
for which ( N e ^/N e ) 2 A 0.5, is 48, 6 and 1 % in the ARs, closed quiet corona, and open 
corona, respectively. To quantify the bimodality of the plasma in the quiet diffuse corona, 
the density added by the hotter component of the N2 model when using the AIA-4 data 
was compared to: a) the electron density computed from the AIA-3 data set, and, b) 
the mean electron temperature computed from the AIA-3 data set. In the first case the 
Pearson correlation coefficient is p = 0.6, while in the later one is p = 0.5, revealing that 
the bimodality is stronger for denser and hotter regions (see Figure IT71). 

Each tomographic cell is threaded by a large number of thermally isolated coronal 
loops and the LDEM obtained at any given voxel is due to the many small-isothermal 
pieces of loops, each one at its own temperature. The fact that the LDEM at any given 
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voxel is bimodal indicates that the typical voxel contains two clearly distinct populations of 
loops: i) ’’warm” loops (described by the cooler component of the N2 model) and, ii) ’’hot” 
loops (described by the hotter component of the N2 model). Those two kinds of loops are 
ubiquitous in the quiet corona, with one or the other being dominant in any given region 
depending on the local average density and temperature. 


Different kinds of physical phenomena have been proposed to explain coronal heating, 
such as dissipation of Alfven waves, nanoflaring, in both ARs (Schmelz et al. 2010; 
Aschwanden & Boerner 2011) and the diffuse quiet corona (Benz & Krucker, 2002; Freij 
et ah, 2014; Uritsky & Davila, 2014; Hahn & Savin, 2014). ft has been suggested that 
nanoflare heating could be mostly operating not directly in the corona but rather at 
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Klimchuk & Brandshaw (2014) have done the same numerically. 


Both studies have found very large discrepancies with the observations (of up to one to two 
orders of magnitude), in terms of line intensities, blue shifts, and blue-red asymmetries. 
They conclude that chromospheric nanoflare can only explain a very small fraction of the 
coronal heating, and that the mechanism should be acting locally in the corona. Moreover, 
based on geometric and energy release considerations, Klimchuk (2015) estimated that a 
typical frequency for nanoflaring in the quiet corona should be of order 200 sec. In such a 
scenario, the stable temperatures of the quiet corona are explained in terms of omnipresent 
high frequency nanoflaring occuring in a constant fashion, with coronal magnetic free 
energy being permanently ’’recharged” as field lines are driven by the constant motion of 
their footpoints. 


Assuming its coronal nature then, nanoflare heating of the homogeneous plasma of 
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the diffuse quiet sun should occur randomly throughout the whole volume of any given 
tomographic voxel, mixing neighboring magnetic threads inside each coronal loop within the 
voxel. In such a scenario, a broad distribution of temperatures is to be expected within each 
magnetic loop (Klimchuk 2006, Aschwanden & Boerner, 2011), which should be detected 
by the LDEM of the tomographic voxel. In fact, since the LDEM provides a time-averaged 
description of the thermodynamics of loops threading the quiet corona, and the result in 
each tomographic voxel is in itself an average of the plasma contained in such a voxel, it 
seems reasonable to expect that the LDEM in each voxel should be very well described by 
broad unimodal distributions, such as the TH or ATH models. The observed bimodality of 
the LDEM seems then difficult to explain in a nanoflare heating scenario. 

The cooler component of the model N2 of the LDEM is an average description of the 
thermodynamics of the warm loops populating the voxel, while the hotter component is an 
average description of the hot loops within the same voxel. In a wave-dissipation dominated 
heating scenario, the reason for having neighboring warm and hot loops within any given 
voxel is to be found in the inhomogeneous chromospheric boundary conditions, where the 
loops are rooted and the injection of energy in waves takes place. Inhomogeneities will lead 
to injection of more energy in waves into some loops (hot ones) and less energy into others 
(warm ones). In this context, we speculate that the observed bimodality of the LDEM may 
be the result of some sort of effective thermodynamical selection effect, allowing only loops 
with very specific characteristics (plus some range of variability) to survive long enough. 

The multimodality is even stronger in ARs, but the tomography is not well suited to 
study such inhomogenoeus and quickly evolving regions. In Section 12.11 we showed how 2D 
DEM analysis can be performed from instantaneous image sets using the same parametric 
inversion technique for the DEM (Equation [8]). We performed such study on a sample AR 
using images in all bands of the AIA-4 set, assuming the N2 parametrization for the DEM. 
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We studied the core of the AR and its surrouding diffuse, selecting regions in a similar 
fashion to a recent study by Narukage et al. (2014), who used a MCMC approach to derive 
DEM models for an AR. Our models in the AR and diffuse region exhibit several similarities 
to Narukage’s results. Firstly, Narukage’s models are typically bimodal, with the values of 
the centroids and amplitudes being similar to our results. Also, both models show a stronger 
bimodality in the core of the AR than in its surroundings, with the centroid temperature 
and area of the second normal component increasing as one gets deeper in the core of the 
AR (see Figure fl8l) . This comparison serves as a validation of the parametric inversion of 
the LDEM that is used in DEMT studies. It also gives support to the bimodal LDEM 
ubiquitously present in the quiet diffuse corona as a real physical property that heating 
models should be able to reproduce, in particular when integrated over the spatiotemporal 
resolution of the tomographic grid. 

F.A.N. acknowledges the CONICET Type I pre-doctoral fellowship that supports his 
participation in this research. The work of E.L. was supported by several NSF and NASA 
grants. 
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Instrument 

-^o lAl 

T 0 [MK] 

FWHM 0 [MK] 

AIA 

171 

0.86 

0.60 


193 

1.50 

0.78 


211 

1.87 

0.87 


335 

2.47 

2.58 

EUVI 

171 

0.98 

0.74 


195 

1.57 

0.82 


284 

2.17 

1.10 


Table 1: Maximum sensitivity temperature (To) and full width half maximum (FWHMo) for 
SDO/AIA and STEREO/EUVI. 
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Fig. 4.— CR-2099. Carrington maps of N e , T m and the score R at 1.075 R Q using data 
of EUVI-B/N1 (left) and AIA-3/N1 (right). In all maps we overplot the boundary between 
magnetically open and closed regions as a thick black curve. ZDAs and AEVs (see text for a 
definition) are indicated as black and white cells, respectively, in the T m maps, and as black 
cells in the N e maps. In the R maps the ZDAs are indicated as black cells. 
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Fig. 5.— Histograms of R at 1.075 Rq in the open region (top panels) and in the closed 
region (bottom panels) using EUVI-B/N1 (left) and AIA-3/N1 (right). 
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Fig. 6.— Comparison of DEMT results in closed magnetic region between EUVI-B/N1 (x- 
axis) and AIA-3/N1 (y-axis) at 1.075 i? 0 . Top panels: scatter plots for N e (left) and T m 
(right), comparing in each panel the results obtained with both instruments. Bottom panels: 
histograms of the ratio of the results shown in the scatter plots. 
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Fig. 7.— Results of inverting the AIA-4 data set using unimodal parametric models. Top 
three panels: Carrington maps of N e , T m and R , at 1.075 R Q . ZDAs and AEVs are marked 
in the same way as in Figure [H Bottom two panels: histograms of the score R separating the 
voxels in the open and closed magnetic region. The three columns correspond to the different 
models. Left: N1 model (3 free parameters), Middle: TH model (4 free parameters), Right: 
ATH model (5 free parameters). See the text for details on the models. 
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Fig. 8.— CR-2099. Carrington maps of N e , T m and R at 1.075 R e using data of AIA- 
4/N2(left) and AIA-4/N4 (right). The boundary between the magnetically open and closed 
region is overploted as a thick black curves in all maps. ZDAs and AEVs are marked in the 
same way as in FigureS] 
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Fig. 9.— Histograms of R at 1.075 Rq in the open region (top panels) and in the closed 
region (bottom panels) using AIA-4/N2 (left) and AIA-4/N4 (right). 
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Fig. 10.— Histograms of the ratio ^ synth )/^( tom ) f or the four bands k of AIA-4, for the 
different parametric models used for the LDEM: N2 (blue), N4 (red), N1 (black), TH (pink), 
and ATH (green). 
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Fig. 11.— Comparition of DEMT results in closed magnetic region between AIA-3 (x-axis) 
and AIA-4 (y-axis). Similar to Figure |6j 
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Fig. 12.— All LDEM models for two sample cells of the magnetically closed corona, one 
belonging to a hotter and denser region than the other. Inversion of AIA-3: N1 (Dahed- 
Blue). Inversion of 4 bands: N2 (solid-blue), N4 (red), ATH (green), TH (pink) and N1 
(black). The values of R , N e and T m for each model is shown in the top region of the Figure. 
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Fig. 13.— Comparison of DEMT results in the closed magnetic region between N2 (x-axis) 
and ATH models (y-axis) for the AIA-4 data set. 
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Fig. 14.— For CR-2099 at 1.075 Rsun we show the relative fraction of plasma (N e j/N e ) 2 
for each component (j = 1, 2) of the N2 LDEM when applied to the AIA-4 data. 
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Fig. 15.— Histograms of ( N ej2 /N e ) 2 at 1.075 R Q , which measures the fraction of plasma 
described by the hotter component of the parametric N2 model. The analysis is shown 
separately for the open region (top panel), the quiet corona in the closed region (middle 
panel), and the active corona in the closed region (bottom panel). 
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Fig. 16.— Top panels: Carrington maps at 1.075 R 0 of the centroids of the first (T 01 ) 
and the second normal distribution (T 0 , 2 )- Center and Bottom panels: Histograms of the 
centroids shown in the top panels in the open and closed magnetic region, respectively. 
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Fig. 17.— Left panel: Scatter plot of N e of AIA-3 versus {N e ^/N e ) 2 at 1.075 R e . Right 
panel: Scatter plot of T m of AIA-3 versus {N e ^/N e ) 2 at 1.075 R Q . In both panels, the blue 
dots corresponds to voxels in the open region and the orange dots correspond to voxels in 
the closed region. In each panel the Pearson correlation coefficient p between both plotted 
quantities is indicated. 
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Fig. 18.— Left panel: EM image based on the parametric inversion of an AIA-4 image 
set of AR NOAA 11089, taken on 2010, July, 25, 12:00 UT. The color boxes, of size 5x5 
pixels, highlight selected subregions in the AR and the diffuse region around it. Right panel: 
Average DEM curve obtained for each color box of the top panel (same color code), using 
the parametric N2 model. 
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